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Abstract. We show that a soliton scattered by an external delta potential splits into 
two solitons and a radiation term. Theoretical analysis gives the amplitudes and phases 
of the reflected and transmitted solitons with errors going to zero as the velocity of the 
incoming soliton tends to infinity. Numerical analysis shows that this asymptotic relation 
is valid for all but very slow solitons. 

We also show that the total transmitted mass, that is the square of the L 2 norm of the 
solution restricted on the transmitted side of the delta potential is in good agreement with 
the quantum transmission rate of the delta potential. 



I. Introduction and statement of results 

A bright soliton is a travelling wave solution, 

(1.1) u(x,t) = Asech(A(x-vt))exp(iip + ivx + i(A 2 -v 2 )t/2), A>0, v GK, 
of the nonlinear Schrodinger equation (NLS): 

(1.2) idtu + \d 2 x u + u\u\ 2 = 

Its remarkable feature is total coherence - see for instance ^T] for a review of theoretical 
and experimental situations in which bright solitons arise. 

Suppose now that we consider a perturbed NLS, that is, the Gross-Pitaevskii equation, 
by adding an external potential: 

id t u + \d 2 u — qS (x)u + u\u\ 2 = 
u(x, 0) = uo(x) 

If as initial data we take a soliton approaching the impurity from the left: 
(1.4) u (x) = e lvx sech(x - x ) , x < , 

then until time to ~ x o/ v , the propagation will still be approximately given by (jl.lj) . Here 
we put A = 1 and <p — 0. Scaling properties of the delta function show that this allows 
general soliton initial conditions. Thus the velocity, v, and the coupling constant, q, are 
the only parameters of the problem. 

For t > xq/v the effects of the delta potential are dramatically visible and as we show in 
this paper they can be understood using the transmission and reflection coefficients of the 
delta potential from standard scattering theory. 
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Figure 1. Numerical simulation of the case q = v = 3, xq = — 10, at times 
t = 0.0,2.7,3.3,4.0. Each frame is a plot of amplitude \u\ versus x. 



For the soliton scattering the natural definition of the transmission rate is given by 



1.5) 



T s q (v) = -\im 



u(t,x)\ 2 dx, / \u(t, x)\ 2 dx = 2 



where on the right we recalled the conservation of the L 2 norm. The reflection coefficient 
is 



;i.e) 



R s q (v) = -\\m 



\u(t, x)\ 2 dx . 



and T s q {v) + R s q {v) = 1. 

The following result is obtained by a numerical analysis of the problem: 

v 2 ( 1 

(1-7) 7» 



v 2 + q 2 



+ o[- 

v 
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In 21 we explain how a weaker rigorous result is obtained in fOJ Theorem 1]. FigEI shows 
the numerical agreement of T q (v) as a function of a = q/v. 

The leading term on the right hand side of (jl.7|) has the following natural interpretation 
in elementary scattering theory, see for instance ^2]- Since we need it below to formulate 
the result about soliton splitting we review the basic concepts. Thus let 

1 d 2 



2dx 2 



q S (x) 



\ 2 /2)u = 0, 
Xx + B±e iXx 



±x > 



H q = - 

and consider a general solution to (H q - 

u(x) = A±e~ 

The matrix 

5(A) : 

is called the scattering matrix and in our simple case it can be easily computed: 

~a\) r q (\y 











i — > 





S(X) 



r q (X) t q {X) 



where t q and r q are the transmission and reflection coefficients: 



'U 



iX 



(I 



iX — q iX — q 

They satisfy two equations, one standard (unitarity) and one due to the special structure 
of the potential: 



1.9) 



\t q (X)\ 2 + \r q 



t q (X) 



The quantum transmission rate at velocity v is given by the square of the absolute value 
of the transmission coefficient (jl.8|) , 

„,2 



1.10) 



T q (v) 



\t q (v)y 



V + q* 

We recall (see [TOI (2.21)]) that if is a smooth function which is zero outside, say, [—a, —b], 
a > b > 0, then 



exp(— itH q )ip(x) \ 2 dx 



T,(v) 



just as in the nonlinear soliton experiment (jl.7|) . 

Hence fll . 7|) shows that in scattering of fast solitons the transmission rate is well approx- 
imated by the quantum transmission rate of the delta potential - see £0 for more on that 
and the comparison with the linear case. 

Our second result shows that the scattered solution is given by a sum of a reflected 
and a transmitted soliton, and of a time decaying (radiating) term. In other words, the 
delta potential splits the incoming soliton into two waves which become single solitons. In 
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Figure 2. The figure illustrates the convergence, as v — * oo, of T^(v) to 
the expected asymptotic value 1/(1 + a 2 ) for a = +0.6, +0.8, . . . , +1.4 (so 
q > 0). It shows that the large velocity asymptotic behaviour in fact takes 
hold by velocity v ~ 3. 



previous works in the physics literature (see for instance [5]) the resulting waves were only 
described as "soliton-like" . More precisely, for t ^> |ar [/v we have 

u{x,t) = u T {x,t)+u R {x,t)+0 L ~ ((t - |x |/f)~ 1/2 ) +0 L 2(v- 2 ), 

u T {x,t) = e^ T e ixv+l{A T- v2)t/2 A T sec\i{A T {x - x - tv)) , 

U R (X, t) = e iVR e -™+i{Al-v*)t/2 AR sech (^ /? ( x + Xq+ ft,)) ; 



1.11^ 



where 



_ , 2\t g (v)\-l, \t q (v)\>l/2 _ 2\r g (v)\-l, \r q (v)\ > 1/2 

0, |t g (v)| < 1/2, 0, |r fl («)|<l/2 



and 



<p T = axgt q (v) + <p (\t q (v)\) + (1- A 2 T )\x \/2v, 
<p R = &rgr q (v) + <p Q (\r q (v)\) + (1 - A^)|x |/2v, 
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Here t q (v) and r g (v) are the transmission and reflection coefficients of the delta-potential 
(see (Oil ). 

The result is illustrated in FigJHJ We can consider A R (q/v) and Aj<(q/v) as nonlinear 
replacements of R q (v ) and T g (v), respectively. Clearly A T +A R ^ 1 except in the asymptotic 
limits q/v — > , oo. Thus if we consider soliton scattering "particle-like" it is nonelastic. 

In FiglHlwe also see the thresholds for the formation of reflected and transmitted solitons: 



Scattering of solitons by delta impurities is a natural model explored extensively in 
the physics literature - see for instance and references given there. The heuristic 

insight that at high velocities "linear scattering" by the external potential should dominate 
the partition of mass is certainly present there. It would be interesting to see if bright 
solitons seen in Bose-Einstein condensates pQ could be "split" using lasers modeled by 
delta impurities 1 . 

In the mathematical literature the dynamics of solitons in the presence of external po- 
tentials has been studied in high velocity or semiclassical limits following the work of Floer 
and Weinstein jS], and Bronski and Jerrard jlj. Roughly speaking, the soliton evolves ac- 
cording to the classical motion of a particle in the external potential. That is similar to 
the phenomena in other settings, such as the motion of the Landau-Ginzburg vortices. 

The possible novelty in (jl.7j) and (jl.ll|) lies in seeing quantum effects of the external 
potential strongly affecting soliton dynamics. The rest of the paper is organized as follows. 
In 321 we present a more detailed discussion of numerical results, and in 313 we outline the 
methods used in our computations. Finally in 3U we discuss weaker but mathematically 
rigorous versions of (jl.7|) and (jl.llj) and give basic ideas behind the proofs. 



We now give numerical evidence for the results presented in 3D the asymptotics (jl.7|) and 
(jl.lip . We stress that the rigorous results of ^U] provide weaker error estimates and hold 
in a limited time range only. We find it very interesting however that some results which 
are theoretically demonstrated, such as the thresholds (J1.12|) or the long time behaviour of 
the the free NLS, are difficult to verify numerically - see § ^2.21 and 12.31 On the other hand 
things which are hard to prove, such as the existence of limits (|1.5|) . seem to be very clear 
in numerical analysis. 

1 That this might be related to the topic of this paper was suggested to the authors by N. Berloff. 



1.12 





2. Numerical results 



6 



J. HOLMER, J. MARZUOLA, AND M. ZWORSKI 




FIGURE 3. A comparison of the linear and nonlinear (soliton) scattering 
rates as functions of a = q/v. 

2.1. Asymptotics for the transmission and trapped rates. We recall from ^ the 
definitions of the transmission and reflection rates 



1 r + oo 

T ? >) = - hm / \u{x,t)\ 2 dx 



1/2 

1/2 

R s g (v) = ^ lim / \u(x,t)\ 2 dx . 



1 

2 t^+oo 



These definitions assume the existence of the limits. The numerical evidence strongly 
supports that the limits indeed exist. We also define the trapped rate 



1 f +1/2 

B s {v) = - lim / \u(x,t)\ 2 dx 

2 <- > +oo < /_ 1/2 



We have 

5>) = 0, g>0. 
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Figure 4. This graph is a log-log plot showing that measured values of 
T*(v) for v = 3, 4, . . . , 10 conform well to the asymptotic formula (J2.3j) . 
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a(a) 


b(a) 




a 


a(a) 


b{a) 


0.6 


0.0415 


2.0786 




-0.6 


0.0441 


2.1076 


0.8 


0.0748 


2.0788 




-0.8 


0.0761 


2.0873 


1.0 


0.1007 


2.0788 




-1.0 


0.1014 


2.0823 


1.2 


0.1147 


2.0778 




-1.2 


0.1151 


2.0798 


1.4 


0.1185 


2.0762 




-1.4 


0.1189 


2.0776 



Table 1. The left table gives the regression coefficients for ()2.3j) for q > 0, 
and the right table for (|2.3j) for q < - see (|2.3j) for the definitions of a(a) 
and b(a), a = q/v. 

However, B*(v) > for q < due to the presence of a bound state at A = —iq for the 
linear operator H q : 

(2.1) <j>(x) = y^q~\e^ . 
The nonlinear problem has a bound state as well 

(2.2) u{x,t) = e lA2 * /2 Asech(A|x|+tanh" 1 (|g|/A)), < A < \q\ . 

This bound state is "left behind" after the interaction. The parameter A depends on the 
initial condition. 
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Figure 5. The first plot shows, for a = —0.6, v = 2, x = —10, the stabi- 
lization of the value of /°q 5 \u(x, t)\ 2 dx (which enters into the definition of 
B s {v)) after the interaction. The second plot shows that the measured value 
of B s q (v) for velocities v = 2.0, 2.5, . . . , 5.0 conforms well to the asymptotic 
formula ()2.4|) . The values of B^(v) that are < e~ 14 cannot be measured with 
adequate precision, and thus the data for \a\ > 1.2 is limited. 
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a 


d(a) 


/(«) 


-0.6 


0.3610 


1.9121 


-0.8 


0.6952 


2.3619 


-1.0 


1.0328 


2.9331 


-1.2 


1.1506 


3.5784 


-1.4 


1.1351 


4.3054 



Table 2. Numerical results for the parameters in (|2.4|) for q < 0. 



In both the q > and g < cases, T|(f) is numerically shown to follow an asymptotic 

(2-3) W~^-4rT> « = -• 

v ; 9 v ; 1 + a 2 v b ^ v 

FigE] is a log-log plot of data for velocities 3 to 10 and a = +0.6 to +1.4. Results of 
linear regression on the subset of the data for velocities 5 to 10 yields the values of a (a) 
and 6(a) reported in Table ^ We see that the value of 6(a) displays little variation with 
a and is approximately 2.07. When the regression is performed on the subset of the data 
for velocities 8 to 10, the values of 6(a) obtained are approximately 2.05 leading us to 
conjecture that the true value of 6(a) is exactly 2, as stated in (J1.7I) . 

The data for q < gives a plot nearly identical to Fig0] owing to the exponential decay 
of the trapped mass B*(v) (discussed below). The values of a(a) and 6(a) obtained from 
the q < data are also reported in Tabled! and we still expect that the true value of 6(a) 
is 2. 

Another feature apparent in Fig0]and Tabled is that the value of a (a) stabilizes asa-> 
oo (note the proximity of the lines for a = 1.2 and a = 1.4 in Fig0]). This feature coincides 
with our analytical result, which establishes a (nonoptimal) bound on the asymptotic of 
f^ 1 / 2 " 1- , independent of a. 

The trapped mass coefficient B*(v), on the other hand, decays exponentially. 
(2.4) B s Jv) ~ d{a)e- f{a)v , a = - 



v 



The second frame of FigE3 which presents data for velocities 2 to 5 and a = —0.6 to —1.4, 
demonstrates that B*(v) conforms well to the formula ([2.4)1 . Linear regression on the data 
yields the values of d(a) and f(a) reported in Tabled 

We see that in constrast to the behaviour of 6(a), /(a) increases with a. This produces 
a numerical road block in studying the asymptotic behaviour further for a = —1.2 and 
a = —1.4. We do not have enough significant digits in our data to measure values of B s q (v) 
less than < e -14 . 
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\qx\ 



For the nonlinear bound state, u(x, t), given by (|2.2|) . we have 

\\u(.,t)\\l 2 = 2(X-\ q \). 

Hence, the behaviour of B*(v) (see (12. 4|) and Table 0J) shows that A approaches q at an 
exponential rate as v — > oo. We note that for A very close to \q\, (|2.2|) is approximately 2 

u(x,t) « e^ 2 */ 2 ^ (l - Ae- A l-I « e^ 2 V2 (l - ^ |g|e" 

which is a multiple of the eigenstate (|2.1j) . 

Given that the interaction with the delta potential is dominated by the linear part of 
the equation, we expect the trapped state will, immediately after interaction, resemble a 
linear eigenstate that will then resolve on a longer time scale to a nonlinear bound state of 
the form ([2.2)1 . It is thus perhaps more appropriate to reverse the heuristics: Given a very 
small amplitude A, if we set 

we obtain that the eigenstate of H q with amplitude A is close to u(x, t) given by (|2.2[) . that 
is, 

(2.5) Ae iq2t/2 e~ lqxl « u(x,t) solving (l2~2l 

It is reasonable to expect that the nonlinear bound state ultimately selected from an im- 
mediate post-interaction eigenstate of the linear operator H q , is "close", in the sense of 
f[2.5|) . to the starting eigenstate of H q , and indeed, the numerics point in this direction. 
The first frame of Fig|S] shows for a typical case (a = —0.6, v — 2) that a stable trapped 
mass is selected within a reasonable amount of time following the interaction and there is 
little evidence of mass being radiated away from the origin. 

2.2. Resolution of outgoing waves. The stabilization of solitons described in (jl.lljl (and 
in a slightly weaker form rigorously in fW\ Theorem 2]) occurs over long time intervals - 
see the comment at the end of £@] 

Hence for the calculation of the amplitudes At and Ar in (jl.llj) we must alter our 
approach. We begin by solving the nonlinear equation (|1.3j) . However, we also measure the 
L 2 difference between the solution at every time step and the expected profile given by 

exp(ivx + it(l — v 2 )/2)[t q (v)sech(x — x — vt) + r q (v)sech(x + x + vt)]. 

Shortly after the time of interaction with the delta potential, we see this difference attains 
a minimum. At this time, we save the computed solution, u(x), and continue to solve 
forward in time. As we solve forward, we compute 

||NLS a (t)u-NLS (t)u||L2. 



2 This is obtained using the two approximations tanhx m 1 — 2e 2x and sechx sa 2e x for x large. 
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Figure 6. A comparison of theory with numerically computed soliton pa- 
rameters An and At- The first plot is for velocity 3; the second, which shows 
a better agreement, is for velocity 10. Note here we are plotting amplitudes, 
as opposed to fraction of the total mass as in FigEl 
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Since this norm remains negligible we switch to the analysis of the simpler solution of the 
(HI , NLS (tK 

More precisely, we truncate the solution on both sides of the delta to give two sets of 
initial data. We then perform a phase shift in order to give each piece a zero velocity. They 
are then embedded at the center of a very large grid with zeroes outside their computed 
ranges. From here, we solve forward on this larger grid using NLSo(t) in order to observe 
the amplitude stabilization predicted by ([l.lljh The grid is chosen large enough so that 
several amplitude oscillations can occur without interference from accumulated errors at 
the boundary. Though the amplitude continues to oscillate, we time average the amplitudes 
until we see stabilization. 

It is these time averages over significantly large intervals that are reported in FigJHl 
We see a clear agreement with (jl.llj) especially for the higher velocity. However, the 
theoretically predicted thresholds for the formation of the reflected and transmitted solitons 
(11.112)) are hard to verify numerically. 

2.3. Confirmation of the free NLS asymptotics for initial data asechx. We now 
turn to the matter of propagating initial data asechx according to id t u + \d 2 u + \u\ 2 u, 
which has been explored analytically via the inverse scattering method in ^UJ Appendix 
B]. FigJ7| reports the results of an experiment with a = 0.8. The first panel depicts the 
time evolution of the amplitude at the spatial origin, \u(0, t)\, and the second panel depicts 
the deviation of the time evolution of the phase at the spatial origin from that of the 
soliton (2a - l)e i(2a ~ 1)2 * /2 sech((2a - l)x). The amplitude appears to be converging to 
the theoretically expected value of 2a — 1 = 0.6 and the phase deviation appears to be 
converging to the expected value of </?(0.8) = 0.045. 

In regard to the phase computation, it should be noted that although this experiment 
was performed on a (x, t) grid of size 15000 x 20000 with spatial extent —600 < x < 600, 
the reported phase deviation is the difference of two numbers on the order of 100 and 
the obtained values are three orders of magnitude smaller. This opens a possibility of an 
inaccuracy of this long-time computation. 

3. Numerical methods 

In this section we outline numerical methods used to produce the results described in §21 
We discretize our equation, 

iu t + u xx + \u\ 2 u — q5o(x)u = 0, 
u(0,x) = u , 

using a finite element scheme in space and the standard midpoint rule in time. Just as 
the equation itself this method is L 2 conservative. A finite difference scheme can also be 
implemented with an approximate delta function, but then convergence must be determined 
in terms of finer meshes as well as more accurate delta function approximations. In finite 
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Figure 7. Two plots depicting the amplitude and phase of u(0,t) for the 
free equation with initial data 0.8sechx. 



element methods the inherent integration allows us to directly incorporate the delta function 
into the discretization of the problem. A similar scheme was implemented without potential 
in j2|, where the blow-up for NLS in several dimensions was analyzed. 
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We review the method here and describe how to discretize the delta potential. Note 
that we require the spatial grid to be large enough to insure negligible interaction with the 
boundary. For the convergence of such methods without potentials see the references in 



We select a symmetric region about the origin, [-R, R] , upon which we place a mesh of 
N elements. The standard hat function basis is used in the Galerkin approximation. We 
allow for a finer grid in a neighbourhood of length 1 centered at the origin to better study 
the effects of the interaction with the delta potential. In terms of the hat basis the problem 
becomes: 

(u t , v) + i{u x , v x )/2 - i(\u\ 2 u, v) + iqu(0)v(0) = 0, 

u(0,x)=u o , u(t,x) = ^Cutyv, 

where (-, •) is the standard L 2 inner product, v is a basis function and u, u are linear 
combinations of the u's. We remark that since u's are continuous the pairing of uv with 
the delta function is justified. 

Since the v's are hat functions, we have created a tridiagonal linear system with one 
contribution to the central element resulting from the delta function. Let h t > be a 
uniform time step, and let 

u n = ^c v {nh t )v , 

V 

be the approximate solution at the nth time step. Implementing the midpoint rule in time, 
the system becomes: 

- Un, v) + iht (((« n+ i + u n )/2) x , v x ) + ih t q(u n+1 (0) + u n (0))v(0)/2 
= ih t (IOn+i +u n )/2\ 2 (u n+1 +u n )/2,v) , m = ^2,a v v, 

v 

By defining 

y n = [u n+ i +u n )/2, 

we have simplified our system to: 

(Vn, V) + i^{(yn)x, V x ) + i^y n (0)v(0) = l^-{\Vn\ 2 Vn, v) + (u n , v) . 

An iteration method from |2] is now used to solve this nonlinear system of equations. To 
wit, 

(y k n + \v) +i^((yt%,v x ) + i^+i (oM o) = i^(\y k n \ 2 yiv) + (u n ,v). 

We take = u n and perform three iterations in order to obtain an approximate solution. 

The measure of success of this method is indicated by the agreement of the computed 
solutions and the exact solutions such as (jl.lj) . and more remarkably by the agreement 
with the inverse scattering method - see ^2.31 above. When the delta potential is present 
(q 7^ 0) the L 2 norm remains essentially constant as predicted by theory. 
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4. Review of theoretical methods 



The theoretical results are proved in jTHj in the case of a positive coupling constant, that 
is for q > 0. The asymptotic formulae (jl.7|) and are obtained rigorously for times 

(4.1) — + v~ 1+e < t < elogv, 

v 

and with the error term (9^2 (t>~ 2 ) replaced by Oip. (t> _ ( 1_3<: )/ 2 ) for any e > 0, provided that 
xq, the center of the initial soliton, satisfies 

xo < —v e . 

The starting point of the argument is the observation that the Strichartz estimates for 
H q hold uniformly in q > 0. Strichartz estimates for dispersive equations [B] describe joint 
space-time decay properties of solutions and are crucial in the control of interaction terms 
in nonlinear equations. 

More precisely, for the problem 

(4.2) id t u(x,t) + \dlu(x,t) - qS (x)u(x,t) = f(x,t) , u(x,0) = u (x) . 
they generalize the well known energy inequality: 

\\U\\ L ^LI < C\\U \\ L 2 + C\\f \\ L l L 2 , 

where, for p ^ oo, 



"Ikfij = I / ( / \u(x,t)\ r dx ) dt 



The Strichartz estimates allow more general exponents p and r and, which is essential to 
us, they hold with constants C = C(p,r,p,f) independent of q: 

\\u\\ LPtL r<C\\u \\ L * + C\\f\\ L p LV 

(4.3) 2 11 2 15 

2<p,r<oo, l<p,r<2, - + - = -, ^ + ^ = 

p r 2 p r 2 

see [TUl Proposition 2.2]. 

We also recall how the reflection and transmission coefficients defined in ^for stationary 
scattering enter in time evolutions: for smooth xjj vanishing outside of [—b, —a], < a < b, 
we have 

r(v)e- itHo [e- lxv ij(-x)](x) + e- itH °[e ixv <ip(x))(x) + e(x,t) , x < , 
t(v)e~ itHo [e ixv ip(x)} (x) + e(x, t) , x > , 

where 

||e(x,t)|| L 2 < -H^VlUa 

v 
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uniformly in t - see [TTH Lemma 2.4]. 

To describe the proof of weaker versions of (jl.7|) and (jl.llj) it will be useful to denote by 

NLS q (t)(p(x) = u(x,t) 

the solution to 

id t u + \d 2 u — qSo(x)u + \u\ 2 u = , u(x, 0) = tp(x) . 

When q > we refer to NLS g (t)</? as the "perturbed nonlinear flow" and when q = 0, as 
the "free nonlinear flow". Similarly exp(—itH q )ip for q > is the "perturbed linear flow", 
and exp(—itH )if is the "free linear flow". 

As discussed in Sect ^ we are interested in 
u(x, t) = NLS q (t)u (x) , u (x) = e lxv sech(x - x ) , v > 1 , x < -v e , < e < 1 . 

The proof of f)1.7|) and (jl.llj) (in the time interval (j4.1j) and with the error term 0^2 (v~ 2 ) 
replaced by (f"^ 1 " 36 ^ 2 )) procceeds in four phases. 

Phase 1 (Pre-interaction). Consider < t < t±, where t\ = \xq\/v — v~ 1+e so that 
xq + vt\ = —v e . The soliton has not yet encountered the delta obstacle and propagates 
according to the free nonlinear flow 

(4.4) u(x, t) = e- ttv2/2 e it/2 e txv sech(x - x - vt) + 0(qe~ vt ), < t < h . 

The analysis here is valid provided v is greater than some absolute threshold (independent 
of q, v, or e). But if we further require that v be sufficiently large so that 

*r 3 /V £ > q/v , 

then 

qe- vi < ir 1 / 2 < v-^' 2 . 
This is the error that arises in the main argument of Phase 2 below. 

Phase 2 (Interaction). Let t 2 = t\+2v~ l+e and consider t\ <t < t 2 . The incident soliton, 
beginning at position — v € , encounters the delta obstacle and splits into a transmitted 
component and a reflected component, which by time t = t 2 , are concentrated at positions 
v e and — v £ , respectively. More precisely, at the conclusion of this phase (at t = t 2 ), 

(4.5) u(x, t 2 ) = t(v)e- it2v2/2 e lt2/2 e lxv sech(x - x - vt 2 ) 

+ r{v)e- it2v2/2 e it2/2 e- ixv sech{x + x + vt 2 ) 

This is the most interesting phase of the argument, which proceeds by using the following 
three observations 

• The perturbed nonlinear flow is approximated by the perturbed linear flow for 
h<t< t 2 . 
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• The perturbed linear flow is split as the sum of a transmitted component and a 
reflected component, each expressed in terms of the free linear flow of soliton-like 
waveforms. 

• The free linear flow is approximated by the free nonlinear flow on t\ < t < t 2 . Thus, 
the soliton-like form of the transmitted and reflected components obtained above is 
preserved. 

The brevity of the time interval [ti,^] is critical to the argument, and validates the ap- 
proximation of linear flows by nonlinear flows. It is here that we used the independence of 
q > in flOp . 

Phase 3 (Post-interaction). Let £3 = £2 + elogf, and consider [£2,^3]- The transmitted 
and reflected waves essentially do not encounter the delta potential and propagate according 
to the free nonlinear flow, 

(4.6) u(x,t) = e- itv2/2 e lt2/2 e lxv NLS {t - t 2 )[t(v)sech(x)](x - x - tv) 

+ e- itv2/2 e lt2/2 e- ixv NLS {t - t 2 )[r(v)sech(x)](x + x + tv) 

+ 0(v-^ 1 - 3e ^), t 2 <t<h 

This is proved by a perturbative argument that enables us to evolve forward a time e log v 
at the expense of enlarging the error by a multiplicative factor of e sv = v e . The error 
thus goes from ti~( 1 ^ e )/ 2 at t = t 2 to u - ( 1-e )/ 2+£ at t = t 3 . 

Phase 4 (Soliton resolution). The last phase uses (I4.fi|) and the following result based 

on inverse scattering method: 

(4.7) 

e *v(«)NLS ((2a - l)sech((2a - 1).)) + L co(H) 1/2 < a < 1 , 
Oi»(r5) 0<a<l/2, 



NLSo(asech) 
where 



f, / sin 2 na \ C , 

v{a) = l log ( 1 + ^hvJ c + 1 /2<«<i, 

see |H3 Appendix B] for the proof and references. The crucial part of the argument involves 
an evaluation of the transmission and reflection coefficients for the Zakharov-Shabat system 
Lip = Xip, 



-iJd x + iJQ, Q = Q(t,x) 



a sechx 
-a sechx 



J 



-1 
1 



That is done by a well known computation [T4*]. |13t Sect. 3. 4] which reappears in many 
scattering theories, from the free S'-matrix in automorphic scattering, to Eckart barriers in 
quantum chemistry. We quote the results: 

r(| + a — iA)r(| - a — i\) . . sin7ra , 
((A) = ^ F(F^ b ^- '^~y r(A) = b(X)t(X) . 
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The pole of t(X) in ImA > determines the soliton appearing in (|4.7|) . We refer to fOJ 
Appendix B] for a detailed discussion. We only mention that the long time asymptotics in 
the case of defocusing NLS [7] show that the error estimate Oloo (t" 1 / 2 ) is optimal and the 
stabilization to the soliton occurs over a long time. 
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